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ABSTRACT 


Calculations  supporting  the  lifting  reentry  communications  systems 
study  are  described  in  this  volume.  The  application  and  interpre¬ 
tation  of  these  calculations  is  prejented  in  Volume  I  of  this  report. 
These  calculations  include  the  study  of  lifting  reentry  trajectories, 
aerodynamic  calculations,  signal  attenuation,  system  margins,  sys 
tern  modifications,  communications  fin  heat  transfer,  magnetic  . 
window,  coolant  injection,  electrophilic  seeding,  and  quasi-optical 
and  optical  systems. 
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i.  INTRODUCTION 


This  report  is  the  second  of  three  volumes  describing  the  analysis  of 
lifting  reentry  communications  systems.  This  volume  contains  the  details 
of  investigations  that  form  the  basis  of  the  argument  in  Volume  I;  Volume  III 
contains  an  extensive  tabulation  of  transmission  and  reflection  coefficients 
fr' -  a  plasma  slab. 

Each  investigation  presented  m  this  report  has  been  specifically  referred 
to  in  Volume  I  as  an  appendix,  Appendices  A  through  J.  These  analyses  are 
intended  to  have  meaning  only  in  the  context  of  the  total  argument  presented 
in  Volume  I,  and  this  relationship  is  emphasized  in  each  section  of  this 
volume . 
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APPENDIX  A 


LIFTING  REENTRY  TRAJECTORIES 


This  appendix  gives  a  short  derivation  of  the  equations  of  motion  of  a 
gliding  reentry  vehicle  and  discusses  the  assumptions  made  to  obtain  the 
trajectories  shown  in  Figs.  1  through  3  of  Volume  I,  Section  II-A-1. 

i.  EQUATIONS  OF  MOTION 

The  flight  geometry  is  illustrated  in  Fig.  A-l.  The  pertinent  symbols 

are 

Cp  =  drag  coefficient 
CT  =  lift  coefficient 

-Li 

D  =  aerodynamic  drag 

g  =  acceleration  of  gravity,  assumed  constant 
H  =  altitude  above  earth  surface 
L  =  aerodynamic  lift 
M  =  Mach  number 
m  =  mass  of  the  R/V 
R  =  earth  radius 

r  =  radial  position  of  vehicle  from  center  of  earth 
Re  =  Reynolds  number 
S  =  aerodynamic  reference  area 
T  =  local  temperature 
U  =  magnitude  of  velocity  vector 
Up  =  orbital  velocity  (r^g/ ^ 

W  -  weight  of  vehicle 


A-l 


a  =  angle  of  attack  of  R/V 
G  -  path  angle 
p  =  local  density 


pg  =  density  at  sea  level 


Fig.  A-l.  Flight  Geometry 


The  definitions 


L  =  SCLpU2/2 


and 


D  =  SCDpU2/2 


lead  to  the  equations  of  motion 


dU 


sin  9  + 


CDS  U2 
^  P  9 


>(>•1 


*U’£3| 


":  ’  .tv1 


=  g  cos  0  - 


U“cos  0  CLSPU 


(A-4) 


The  initial  conditions  are  t  =  0,  U  =  Uq  ,  0  =  0Q,  and  H  =  HQ. 

The  relation  p(H)  may  be  taken  from  a  standard  atmosphere  table  or  be 
approximated  by  the  isothermal  atmosphere 


=  exp[-(g/RT)H] 
ps 


(A-5) 


The  dependence  of  and  on  a,  Re,  and  M  must  be  determined  for  the 
particular  geometry  under  consideration.  The  altitude  of  the  R/V  as  a  func¬ 
tion  of  time  is  given  by 


H  ’  H»  •  f 


U  sin  0  dt 


(A-6) 


2.  EQUILIBRIUM  TRAJECTORY 

An  equilibrium  lifting  reentry  trajectory  is  characterized  by  a  small 
path  angle  and  a  small  rate  of  change  in  the  path  angle.  Under  these  circum¬ 
stances,  the  equations  of  motion  become 


-cds  yf 

m  ^  2 


(A-7) 


u  r 

U0  ~  [i  +  r(C 


^S/m)(p/2) 


(A-8) 


If  we  introduce  p  -  p  exp(  -(g/RT  t Hj ,  the  velocity  can  be  computed  as  a  function 

8 

of  altitude  from  Eq.  (A-8).  Integration  of  Eq.  (A-7)  will  give  t(U).  For 
r  =  Tq  =  constant,  the  integration  yields 


t 


t  3 

i  2g 


L 

D 


i  11 


[1  +  (U./U0)]  ll  -  (U/UQ)] 

11  -  (Uj/Uq)]  UTTuTuJI 


(A-9) 


Trajectory  1  of  Fig.  2,  Volume  I,  has  been  computed  from  Eqs.  (A-8) 
and  (A  9)  using  the  ARDC  Atmosphere  Tables.  1  Trajectory  1  corresponds 
to  an  entry  withUQ=  25.  6  kft/sec,  0=0,  and  a  ballistic  parameter  W/CLS  of 
200  ib/ft2. 

3.  NONEQUILIBRIUM  TRAJECTORY 

As  illustrated  by  Fig.  1,  Volume  I,  a  nonzero  initial  path  angle  leads  to 
oscillations  of  the  trajectory  in  the  altitude  vs  velocity  plane.  To  consider 
the  increased  depth  of  blackout  encountered  during  the  low-altitude  high- 
velocity  pullout,  we  have  defined  an  idealized  and  rather  extreme  nonequilib¬ 
rium  trajectory  which  is  the  envelope  of  trajectories  computed  by  numerically 
integrating  Eqs.  (A-3)  and  (A-4)  for  8q  *  0  and  W/C^A  =  200  lb/ft  .  The 
corresponding  time  scale  was  determined  by  observing  that  the  time  elapsed 
from  the  reentry  point  is  approximately  independent  of  altitude. 


^  The  numerical  integration  of  Eos.  (A-7)  and  (A-8)  was  originally  performed 
by  the  General  Dynamics  Corp.  ,  Fort  Worth,  and  the  results  made  available 
to  us  by  Mr.  W.  C.  Melton  of  ESTO,  Aerospace  Corp. 


APPENDIX  B 


AERODYNAMIC  CALCULATIONS 


This  appendix  is  a  compilation  cf  results  obtained  by  the  methods 
described  in  Volume  I,  Section  II-A-3. 

1.  OUTER  IN  VISCID  FLOW 

Figures  B-l  through  B-4  show  the  plasma  frequency,  degree  of  ioniza¬ 
tion,  and  collision  frequency  in  equilibrium  air.  The  plasma  frequency  and 
degre.  of  ionization  are  taken  from  work  by  Bleviss  (Ref.  B-l),  which  was 
based  on  the  results  of  Logan  and  Treanor  (Ref.  B-2),  The  collision  frequency 
was  obtained  from  Bachynski,  et  al.  (Ref.  B-3). 

Figures  B-5  through  B-8  show'  inviscid  shock  layer  properties  for  cones, 
wedges,  and  axis^mmetric  stagnation  points.  These  figures  were  adapted  as 
follows:  Fig.  B-5  'rom  Ref.  B-4;  Fig.  B-6  from  Ref.  B-5;  Fig.  B-7  from 
Ref.  B-6;  and  Fig.  B-8  from  Refs.  B-7  and  B-8. 

2.  VISCOUS  BOUNDARY  LAYER 

The  boundary  layer  properties  were  obtained  from  an  existing  computer 

program  that  is  described  in  Ref.  B-8.  The  maximum  plasma  frequency  in 

the  boundary  layer  on  cones  and  wedges  is  presented  in  Figs.  B-9  through 

B-12.  Figures  B-13  through  B-16  show  the  collision  frequency.  It  is  to  be 

noted  that  these  properties  differ  by  as  much  as  50°,o  from  the  best  available 

data  (see  Volume  I,  Ii-A-3),  The  normalized  boundary  layer  thickness  6/(x)- 

for  wedges  is  shown  in  Figs.  B-17  through  B-20  (The  boundary  layer  thickness 

in  inches  is  6,  and  the  distance  along  the  surface  of  the  vehicle  measured  from 

the  nose  in  feet  is  x.  );  the  conical  boundary  layer  thickness  is  less  by  a  factor 
-1/2 

of  (3)  .  Figure  B-21  presents  representative  values  of  the  Reynolds  number 

at  the  boundary  layer  edge. 
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Fig.  B-5:  Plasma  Frequency  and  Electron  Collision 
Frequency  in  Equilibrium  Air  Behind  Normal  Shocks 
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Fig.  B-io.  Collision  Frequency 
t  Peak  Boundary  Layer  Condition: 
Shock  Angle  -  20  deg 


Fig.  B-15.  Collision  Frequency 
at  Peak  Boundary  Layer  Condition: 
Shock  Angle  =  30  deg 
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Fig.  B-19.  Wedge  Boundary 
Layer  Thickness  (Laminar): 
Shock  Angle  =  30  deg 
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Fig.  B-20,  Wedge  Boundary 
Layer  Thickness  (Laminar): 
Shock  Angie  -  20  deg 
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APPENDIX  C 

SIGNAL  ATTENUATION 

Signal  attenuation  calculated  as  indicated  in  Volume  I,  Section  II-C-I,  } 

for  cones,  wedges,  and  stagnation  points  are  presented  in  Tables  C-I  through  \ 

f 

C-X1II.  The  calculations  were  made  for  five  signal  frequencies  (0.  25,  0.  5,.  f 

t.  0,  5.  0,  10  GHz)  for  representative  points  along  both  equilibrium  glide  and  J 

transient  trajectories  (see  Volume  I,  Section  II-A-i).  Results  are  presented  I 

I 

according  to  the  following  tabulation:  J 


Configuration 

Body  Station 

Flow  Field 

Table 

Stag.  Pt. 

nose  (3-in.  radius) 

Shock  Layer 

C-I 

Sharp  Wedge 

(1  ft 

Shock  Layer 

C-II 

(50") 

15  ft 

Sharp  Wedge 

(i  ft 

(Shock  Layer 

C-lII,  V 

(30°,  40°) 

(5  ft 

(Boundary  Layer 

C-IV,  VI 

Sharp  Wedge 

it  ft 

Boundary  Layer 

C-VII 

(20°) 

(5  ft 

Sharp  Cone 

(i  ft 

Shock  Layer 

C-VIII 

(50") 

(5  ft 

Sharp  Cone 

(1  ft 

(Shock  Layer 

C-IX,  XI 

(30°,  40c ) 

(5  ft 

(Boundary  Layer 

C-X,  XII 

Sharp  Cone 

it  ft 

Boundary  Layer 

C-XIII 

(20°) 

(5  ft 

Table  C-XIV  presents  the  characteristic  incubation  distance  for  electrons 
behind  plane  shock  waves  in  air  (Ref.  C-i)  compared  with  the  appropriate 
shock  layer  thicknesses.  These  data  have  been  used  in  the  discussion  of  the 
importance  of  nonequilibrium  effects  in  Volume  I,  Section  U-C-2. 


1 


C-I. 


C-II.  Shock  Lay 


Table  C-I1I.  .Shock  Layer  Attenuation,  40-deg  Wedge 


Table  C- V.  iihock  Layer  Attenuation,  30-deg  Wedge 


Table  C-VI.  Boundary  Layer  Attenuation,  30-deg  Wedge 


Ml  I  llx  I  t. 14 


Table  C-V1I.  Boundary  Layer  Attenuation,  20-deg  Wedge 
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Table  C-VIII.  Shock  Layer  Attenuation,  50-dey  Cone 


'  ►  i  !  :  I  ;  u  *  i  s  « 


5  S  S  3  S  2 


a  s  3  s  S 


:  3  ;  :  ;  ;  j  i  t 


A  <i  M  4  f 

*  A  A*  •• 


r  * 


s  :  :  :  j  s  i  s  :  :  j  s  -  isi  'i:s 


N  ♦  4  f 


»  -t  4  «  4 


5  S  4  r  ?  t  5 


J«*i«AKeoc>w^Np» 
r  «  »  o  o  -  •  «  f.  > 


t  ;  i  ;  d  : 


o  o  o  o  o 


oooo«oooooooo 


4  -  a  «  *  v  4 

«  M  «  «• 


£  C  S  * 


WOO  090000 

• 

;i  »««»»»»  ’ 


3  u  u  3  :  ;  :  :  : 


dbeooooeboooo 


VVXXKXKXXM  K 


#  -  •»  r  •  ► 


*  x  X  X  X  i:  X 
m  o  •  •  «ft  •  • 
t  *  «  •  «  4  • 


X  u 

f  •  «  e  e  o 

tf£  .;  «  «  « 

-S  M  K  H  »* 


d  fl  d  d  ♦  4  4 


S  S  3 

<*4  r*i  n 


S*  “  8  2 


n  U  S  H  !  1  U  H 


°asss2?sasss 


5  5  l  l  5  3  5 


Table  C-IX.  Shock  Layer  Attenuation,  40 -deg  Cone 


Table  C-X.  Boundary  Layer  Attenuation,  40-dey  Cone 
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Table  C-XI.  Shock  Layer  Attenuation,  30-deg  Cone 


m  I  14.0  I  o.ot  j  u»ir  |  mi  |  s.ss  |  o.oo  |  •. »  I  o.o  I  o.o  I  o.o  I  o.o  1  o.o  I  o.o  I  o.o  1  0.0 


Table  C-XII.  Boundary  Layer  Attenuation,  30>deg  Cone 
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Table  C-XUI.  Boundary  Layer  Attenuation,  20 -deg  Cone 
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Table  C-XIV.  Shock  Layer  Thickness  vs  Ionization  Distance  for  Cones  and  Wedg 
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APPENDIX  D 


SYSTEM  MARGINS 


As  discussed  in.  Volume  I,  Section  II-D,  the  free  space  system  margin, 
which  is  defined  as  the  amount  of  plasma  attenuation  that  the  reentry  vehicle 
system  can  tolerate  while  it  effectively  transmits  information,  is  a  function 
of  frequency-  The  calculations  which  lead  to  the  frequency  dependence 
exhibited  in  Fig.  23  of  Volume  I  are  presented  here. 

The  free  space  system  margin  M  is  given  by 

M  =  ■'VlC”S.' 'N>0/(S/N)T'j  (D-l) 

where  (S/N)Q  is  the  signal-to-noise  ratio  neglecting  plasma  attenuation  and 
(S/N)t  is  the  threshold  signal-to-noise  ratio.  If  we  take  the  noise  figure  of 
the  receiver  to  be  unity,  the  link  equation  may  be  written 


(D-2) 


where  \q  is  the  free  space  wavelength,  R  is  the  range,  P^,  is  the  transmitter 
power,  G-j,  is  the  transmitting  antenna  gain,  Gj^  is  the  receiving  antenna  gain, 
Af  is  the  rf  bandwidth,  L  represents  miscellaneous  losses,  k  is  Boltzmann's 
constant,  and  Tg  is  an  effective  noise  temperature  given  by 

Te  =  290*K  +  TA  (D-3) 

where  is  the  effective  antenna  temperature. 


D-l 


As  discussed  in  Volume  I,  parameters  were  chosen  as  follows: 


(S/N)t  =  10 

R  =  100  mi 
Xq  =  0.  03  to  300  m 
PT  =  l  W 


f  =  300  kHz 
L  =  10 

The  receiving  antenna  gain  and  effective  antenna  temperature  are 
dependent  on  frequency,  and  this  dependence  is  illustrated  in  Figs.  D-l  and 
D-2.  These  graphs  are  intended  to  be  representative  rather  than  to  provide 
ultimate  design  criteria.  The  data  points  of  Fig.  D-l  were  obtained  from 
manufacturer's  catalogues  and  published  characteristics.  The  data  of  Fig.  D-2 
w«re  otiained  from  Refs.  D-l,  D-2,  and  D-3.  In  the  frequency  range  below 
10  MHz  the  noise  is  largely  atmospheric,  whereas  in  the  range  between  10  MHz 
and  500  MHz  the  noise  is  largely  galactic. 

The  result/  of  the  system  margin  calculation  are  given  in  Fig.  D-3.  As 
discussed  in  Volume  I,  system  margins  for  other  transmitter  power  levels, 
receiver  bandpasses,  antenna  gains,  and  ranges  cat  easily  be  obtained  from 
Fig.  D-3  and  Eq.  (D-2). 
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Fig.  D-l.  Antenna  Cain 
vs  Frequency 
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APPENDIX  E 

SYSTEM  MODIFICATION 

Calculations  pertinent  to  the  evaluation  of  the  feasibility  of  increasing 
the  rf  power  of  lifting  reentry  communications  systems  are  presented  in  this 
appendix.  The  application  of  the  results  is  presented  in  Sections  III-B, 

III-D,  and  III-E  of  Volume  I. 

1.  POWER  AMPLIFIER  WEIGHT  PENALTY 

The  weight  of  a  power  amplifier  over  the  range  of  conventional  telemetry 
frequencies  may  be  considered  to  be  directly  proportional  to  both  the  power 
output  and  the  operating  frequency  (Ref.  E-l).  At  250  MHz, a  6-oz  water- 
cooled  tube  (EIMAC  4W300B)  delivers  about  200  W  of  rf  power  with  an  input 
of  about  300  W  of  dc  power.  A  power  supply  capable  of  delivering  300  W  for 
a  10 -min  period  should  weigh  about  3  lb  (this  estimate  includes  batteries  and 
converter).  The  water  required  to  dissipate  100  W  for  10  min  is  less  than 
0.  1  lb.  Al'owing  1  lb  for  the  tube  and  associated  components  and  3  lb  for  the 
power  supply,  we  estimate  the  total  weight  for  the  200-W  amplifier  to  be 
4  lb.  Using  this  estimate  for  normalization,  we  derive  the  relation 

W  =  PQf(l.  25  X  104)'1  lb  (E-l) 

where  W  is  the  weight  penalty,  PQ  is  the  output  power  in  W,  and  £  is  the 
frequency  in  MHz. 

We  observe  that  an  EIMAC  X-1134  Telemetry  Amplifier  Package  delivers 
20  W  of  CW  rf  power  at  8000  MHz.  The  package  weight  is  10  lb,  indicating  a 
total  weight  of  11  lb  if  we  allow  1  lb  for  a  200-W  28-V  source.  Our  weight- 
penalty  formula  gives  13  lb  for  such  a  system,  and  this  is  taken  to  indicate 
that  we  may  use  the  formula  to  make  sufficiently  accurate  estimates. 


E-l 


2. 


BREAKDOWN  LIMITATIONS  FOR  PULSED  SOUR CES 


Radio  frequency  breakdown  limitations  for  CW  sources  have  been 
discussed  in  Volume  I,  Section  II-B-4.  Here  we  determine  whether  break¬ 
down  limitations  can  be  reduced  by  using  a  source  that  is  pulsed  in  such  a 
way  that  average  power  remains  constant  while  the  duty  cycle  is  reduced. 

For  pulsed  rf  energy,  breakdown  occurs  when  the  ionization  rate 
multiplied  by  the  pulse-length  r  exceeds  a  certain  constant  (Ref.  E-2). 

Since  the  pulse-length  is  inversely  proportional  to  the  peak  power,  the  power 
is  proportional  to  the  square  of  the  electric  field  E,  and  the  ionization  rate  is 
a  function  of  the  electric  field,  we  conclude  that  the  condition  for  avoiding 
breakdown  may  be  written 


v,(E) 

R  = - y—  <  constant  (E-2) 

el 

Table  E-I,  which  is  based  on  calculations  presented  in  Ref.  E-2,  indicates 
that  the  quantity  R  defined  by  Eq.  (E-2)  is  an  increasing  function  of  peak 
power  and,  therefore,  a  decreasing  function  of  pulse -length.  We  conclude 
that  breakdown  limitations  become  more  severe  for  a  given  pulse  of  electro¬ 
magnetic  energy  as  the  pulse-length  is  reduced  (see  Volume  I,  Section  III-B). 


Table  E-I.  Behavior  of  Breakdown  Parameter  R 


E  (relative) 

(E)  (relative) 

2 

E  (relative) 

R  (relative) 

2 

1.  5 

4 

0.  4 

3 

15 

9 

1.  7 

4 

40 

16 

2.  5 

5 

150 

25 

6 

10 

3000 

100 

30 
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APPENDIX  F 


COMMUNICATIONS  FIN  HEAT  TRANSFER 


In  Fig.  25  of  Volume  I,  Section  III-C,  results  of  heat  transfer  calcula¬ 
tions  are  presented  to  assess  the  feasibility  of  the  communications  fin  con¬ 
cept.  The  relationships  used  in  these  calculations  are  reproduced  here.  The 
notation  is  defined  in  the  referenced  literature  (Ref.  F-l). 


1.  STAGNATION-LINE  HEAT  TRANSFER 

The  heat  transfer  to  the  0.  25 -in.  radius  leading  edge  of  the  communica¬ 
tion  fin  has  been  calculated  for  the  entire  nonequiiibrium  trajectory  using  the 
following  relationship: 


0.  576 


0.  44, 


lst 


Pr 


^<Ve>  <*w»w> 
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0.  5 


{■ 


<hr  •  Vl,  +  <Le°'52  *  *>fC 


LD] 

e  J 


(F-l) 


This  expression  is  a  modification,  by  Kemp,  Rose,  and  Detra  (Ref.  F-2), 
of  the  original  relationship  given  by  Fay  and  Riddell  (Ref.  F-3).  Conditions 
at  the  wall  were  arbitrarily  chosen  to  correspond  to  Tw  =  273?K. 

2.  FLAT-PLATE  HEAT  TRANSFER 

The  heat  transfer  to  the  side  of  the  fin  and  1  ft  from  the  leading  edge  has 
been  calculated  using  the  relationship 


1.  383^(Re*)1/2(Pr*)1/3(h.  -  h  J 


r  w 


(F-2) 


F-l 


where  k*,  Re*,  and  Pr*  are  evaluated  at  a  temperature  T*  given  by 


pideal 


0.  28h  +  0.  5h  +  0.  22h 
_ e  w _ r 

0.  24 


(F-3) 


The  ideal  gas  equation  of  state  and  the  Sutherland  viscosity  law  are  assumed. 
Wall  conditions  were  chosen  to  correspond  to  Tw  =  1300*K.  Heat  transfer 
rates  calculated  were  compared  with  the  calculations  of  Miles  and  Waldman 
(Ref.  F-4)  and  were  found  to  be  in  good  agreement. 
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APPENDIX  G 


MAGNETIC  WINDOW 


Formulas  and  sample  calculations  pertinent  to  magnetic  window 
systems  are  presented  here.  The  application  of  these  results  to  the  evaluation 
of  magnetic  window  as  a  reentry  communications  system  is  presented  in 
Volume  I,  Section  III-D. 

1.  DERIVATION  OF  FORMULAS 

a.  COIL  GEOMETRY 

We  envision  a  coil  composed  of  copper  or  aluminum  wire  wrapped 
around  a  circular  antenna  aperture.  Figure  G-l  is  a  cross-sectional  view 
of  the  coil  showing  the  average  diameter  d,  the  cross-sectional  width  a,  the 
aperture  diameter  b,  and  the  distance  to  the  sheath  y. 


CURRENT  ASSUiflED 
CONCENTRATED  HERE 
IN  CALCULATION 


Fig.  G-l.  Coil  Geometry 

b.  CURRENT  REQUIREMENT 

The  total  current  IT  in  a  cross  section  of  the  coil  is  equal  to  the 
number  of  turns  in  the  coil  multiplied  by  the  current  at  the  input  terminals. 

The  total  current  requirement  is  linearly  related  to  the  magnetic  fie.^  strength 
required  in  the  sheath.  For  our  purposes,  this  relationship  is  approximated 


by  assuming  that  the  equivalent  uniform  magnetic  field  of  the  coil  has  the 
same  dependence  on  the  displacement  y  as  the  axial  field.  With  this  approxi¬ 
mation,  we  obtain 


3-1) 


and  to  complete  the  specification  of  the  current  requirement  we  estimate 


d 


(G-2) 


In  the  above  expressions,  B  is  the  average  field  in  the  sheath,  Bq  is  the 
average  field  in  the  plane  of  the  coil,  and  is  the  magnetic  permeability  of 
free  space.  (MKS  units  are  used  for  all  quantities  except  weights  which  are 
given  in  pounds.  ) 

c.  SYSTEM  WEIGHT  AND  APERTURE  DIAMETER 

The  total  system  consists  of  the  coil,  the  batteries,  and  the  cooling 
system.  If  we  assign  a  weight  to  the  coil- supporting  structure  equal  to  10% 
of  the  coil  weight,  and  a  weight  to  the  cooling  system  equal  to  twice  the  weight 
of  the  water  required  to  dissipate  ohmic  heat  losses,  we  obtain 


W 


=  nd 


1.  IpA  + 


(G-3) 


where 

W  =  system  weight,  lb 

3 

p  =  specific  density  of  conductor,  lb/m 
A  =  =  cross-sectional  area  of  coil,  (see  Fig.  G-i) 


0-2 


r  =  resistivity  of  coil  material,  ehm-m 

t  =  total  time  of  operation,  sec 

L  -  specific  he?!  of  vaporization  of  water,  J/lb 

B  =  specific  energy  capacity  of  batteries,  J/lb 

The  value  of  cross -sectional  area  A  is  taken  as  the  value  that  minimizes  the 
system  weight.  The  minimization  procedure  yields 

W  =  2.  ?.  x  10"2  lTd(rpt)i/2  lb 

A  =  3.  2  x  10'3  IT(rt/p)i/2  m2 

b  =  d  -  A1/2  m  (G-4) 

We  note  that  system  weight  is  direetly  proportional  to  the  magnetic  field  and 
the  square  root  of  the  on-time. 

The  voltage  input  V.ft  and  the  current  input  lin  are  determined  by  the 
number  of  turns  N  in  the  coil;  the  number  of  turns  is  determined  by  the  size 
of  the  coil  wire.  If  the  cross -sectional  area  of  the  wire  is  <r,  we  have 

N  =  A/<r 

V.n  =  NIx(rnd/A) 

I*n  =  I,,/N  {G-5) 
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SAMPLE  CALCULATIONS 


Representative  requirements  for  a  magnetic  window  system  capable 
of  limiting  signal  attenuation  to  40  dB  at  a  station  l  ft  from  the  apex  of  a 
wedge  with  a  50-deg  angle  of  attack  have  been  presented  in  Volume  I, 

Section  III-C-3-b.  The  estimated  parameters  of  this  system  were 

d  =  4  in. 

y  =  i  in. 

B  -  0.  16  Wb/m2  (ItOQ  G) 
t  =  500  sec 

We  took  the  following  parameters  as  representative: 

L  =  106  J/lb 

B  «  1.  i  X  105  J/lb  (30  W-h/lb) 
o 

r  =  2.  2  x  10  •  ohm-m  (copper) 

Q 

r  =  3*  9  x  i0  ohm-m  (aluminum) 

p  =  2  X  104  lb/m"*  (copper) 

3  3 

p  =  6  x  10  lb/m  (aluminum) 

=  4ir  X  to'7  H/m 

Using  these  values  in  Eq.  (G-4),  we  find  that  for  a  copper  coil  the  system 
weight  is  about  17  lb, and  the  aperture  diameter  is  about  2.  5  in.  For  an 
aluminum  coil,  the  system  weight  is  about  12.  5  lb,and  the  aperture  diameter 
is  about  1.  5  in. 

As  indicated  in  Volume  I,  Section  IV-B,  the  evaluation  of  the  magnetic 
field  approximations  and  aperture  effects  is  designated  as  a  task  for  a  future 
research  effort. 


APPENDIX  H 


COOLANT  INJECTION 


The  various  calculations  pertinent  to  the  evaluation  of  coolant  injection 
as  an  alleviation  technique  as  described  in  Volume  I,  Section  Ill-E,are 
collected  here. 

1.  PROPERTIES  OF  COOLANTS 

Figure  H-l  is  adapted  from  Ref.  H-l.  The  data  in  the  figure  are 
discussed  in  the  argument  for  the  selection  of  water  as  an  inject«*nt  in 
Volume  I,  Section  III-E-1. 


Fig.  H-l.  Cooling  Effectiveness  as  a  Function 
of  Temperature  for  Various  Coolants 


H-l 


As  pointed  out  In  Volume  I,  Section  III-E-2-b,  we  use  a  technique 
identical  to  that  of  Ref.  H-l. 

Consider  injection  into  a  streamtube  as  indicated  in  Fig.  H-2. 


Fig.  H-2.  Ouasi  One-Dimpnsional  Mixing  Process 

The  properties  are  assumed  uniform  across  the  cross  sections  of  the  stream- 
tube  at  Stations  1  and  2.  The  conservation  equations  are  (we  assume  551  0) 

PlUlAl  =  *2U2A2  (H-” 

piUlAl  +  *”c  =  p2U2A2  (h-2) 

P1U1AI  [hl  +HUl)]+  *chc  *  p2°2A2[h2  +  r(Ui)]  ,H'3) 

where  p  is  the  density;  U  is  the  velocity,  h  is  the  enthalpy,  mc  is  the 
injection  rate,  and  A  is  the  cross-sectional  area.  V/e  now  specifically  define 


H-2 


c 


=  mass  fraction  of  coolant  at  Station  2 


p  =  density  of  coolant  at  Station  Z 
c2 


p  =  density  of  air  at  Station  2 
2 

mc 

Cq  =  —  -g  =  mass  ratio  of  injectant  to  incoming  flow 
Noting  that  C  =  Cq/(1  +  Cq),  we  obtain 


h2  =  [1/(1  +  CQ)]  [ht  +  C  (Uj  /2)  4  CQhc 


(H-4) 


Thus,  for  a  given  state  at  Station  1  and  an  injection  ratio  Cq,  we  may 
calculate  the  enthalpy  at  Station  2,  and  the  pertinent  properties  of  the  mix¬ 
ture  may  be  deduced.  The  electron  density  is  obtained  from  data  for 
equilibrium  air  as  a  function  of  the  temperature  and  the  air  density.  The 
collision  frequency  is  obtained  (approximately)  from  the  data  for  equilib¬ 
rium  air  (Appendix  B,  Fig.  B-4),  where  it  is  assumed  that  the  electron- 
neutral  cross  section  for  helium  is  seven  times  greater  than  that  of  air  on 
a  unit  mass  basis.  Hence,  the  ordinate  in  Fig.  B-4  is  interpreted  as 


v(po/peff,#  where 


eff 


P 

~P 


?  •  4) 


I  +  6C 

r-  c 


(H-5) 


and  the  parameter  p/pQ  is  interpreted  as  p  fpQ.  This  is  an  adequate  approxi¬ 
mation  at  the  low  electron  densities  of  interest,  since  electron-ion  Interactions 
do  not  contribute  significantly  to  the  total  electron  collision  frequency. 

The  calculation  of  system  weight  requirements  is  described  in  Volume  I, 
Section  III-E-2 . 


3. 


ELECT RQNEGATIVES  AS  CATALYSTS 


The  relaxation  time  of  air  in  the  presence  of  an  electronegative  gas  is 

given  by  (f^c^n where  <r^  is  the  attachment  cross  section,  is  the 

mean  thermal  speed  of  the  electrons,  and  n-,  is  the  number  density  of  the 

-lb  2 

electronegative  particles.  For  SF^,  is  of  the  order  of  10  cm 
(see  Appendix  I).  For  a  typical  reentry  case  (say,  a  local  velocity  about 
15  kft/sec  and  a  local  pressure  of  about  0.01  atm)  the  equilibration  time  should 

be  about  50  psec.  To  obtain  a  relaxation  time  of  5  psec  n_  should  be  about 

14  ^ 

10  /cc.  Since  SF/  is  about  ten  times  heavier  than  water,  the  water  weight 

o  .17 

is  about  equal  to  the  weight  of  air,  and  the  air  density  is  about  10  / cc,  we 

conclude  that  the  injected  mixture  should  contain  about  a  1%  by  weight 

concentration  of  SF^  (see  Volume  I,  Section  III-F-2-b). 

4.  ESTIMATE  OF  DROPLET  SIZE 

The  following  empirical  expressions  for  droplet  sizes  in  breakup  of 
liquid  jets  in  air  streams  are  taken  from  Ref.  H-2: 


n/Do  =  3.  9(We/Re) 


0.  25 


D  /D  =  22.  3( We / Re ) 
max  o 


0.  29 


(H-6) 


(H-7) 


where 

D  =  volume  average  droplet  diameter 

.0  =  maximum  droplet  diameter 

max 

D  -  orifice  diameter 
o 

W.  = 

Re  =  D  V  /v 
o  s 

<r  ~  surface  tension  of  liquid 


H-4 


p  -  ai»--stream  density 
8 

V  -  air  jtream  velocity 
8 

v  =  liquid  kinematic  velocity 
For  water 

<r  =  65  dyn/cm 

v  =  4.  7  x  10"^  g /sec -cm  (at  60  *C) 


D  =  0.  01  in. 
o 

5 

V  :  4.5  X  10  cm/sec 

s 

p  =  1. 29  x  10"5  g/cc 
we  obtain 

D  =  4.  45  p 


D  =  10.  8  p 
max 


As  discussed  in  Volume  I,  Section  III-E-2-b,  we  conclude  that  a  reasonable 
estimate  for  the  droplet  size  in  the  air  stream  is  from  5  to  10 p. 

5.  EVAPORATION  OF  WATER  DROPLETS 

Free  Molecular  Flow 

For  i  droplet  of  water  in  air 


,2  1  .  d  ,  ,3. 

qd  =  -  (d  ) 


(H-8) 


where  q  is  the  energy  flux,  d  is  the  diameter  of  the  droplet,  p  is  the  specific 
density  of  water,  and  L  is  the  energy  per  unit  mass  required  for  evaporation 
of  water.  We  readily  obtain 

^d qI-'P  I  j-x 

t  =  — 2—  (1  -  f1'/  (H-9) 

e  q 


H-5 


where  tg  is  the  evaporation  time,  dQ  is  the  initial  drop  diameter,  aid  f  is 
the  ratio  of  the  unevaporated  water  weight  to  the  initial  water  weight.  The 
energy  flux  q  is  given  by 

q  =  Z  n.c.E.  (H-10) 

i 

I 

and 


£.  =  2kT  +  -  (H-l  1 ) 

where  the  summation  is  over  species,  n.  is  the  particle  number  density, 
c.  is  the  thermal  velocity,  E.  is  the  energy  per  particle,  and  <f> ^  is  the 
energy  obtained  from  recombination  processes.  1  ak.ing  5.  1  eV  as  the 
recombination  energy  for  oxygen,  9.  8  eV  as  the  recombination  energy  for 
nitrogen,  and  the  concentrations  listed  in  Table  H-I,  we  obtain  the  results 
listed  in  Table  H -II. 


Table  H-I.  Species  Concentration  in  Air 


T.  *K 

[Nj] 

[N] 

[o] 

5950 

0.  34 

0.  37 

0.  28 

4000 

0.  65 

- 

0.  35 

Table  H-H.  Evaporation  Times,  Free  Molecular  Flow 


Evaporation,  % 

te,  psec 

5950°K  rate 

4CO0*K  rate 

70 

0.  062 

0.  27 

80 

0.  078 

0.  34 

90 

0.  100 

0.  44 

100 

0.  190 

0.  55 

H-6 


b.  Continuum  Flow 

For  the  continuum  case,  the  energy  /lux  is  given  by  {see  Ref.  H-; 
references  cited  therein) 


N-.  k,(h  -h  ) 
Nu  f  m  w 


p,  w 


where 


ZC  JL 

N  .  =  — -  - 

Nu  C-  -  7, 


p,  v,  f(h  -  h 
r  m  w 


(C  r  h  -  h  \ 
p,m  / 


and 


NNu 

kf 

hm 

h 

w 

C 

p.  m 
'p.  v.  f 


Nusselt  number 

thermal  conductivity  of  air  at  Tf  =  |  <Tdrop  +  Tmlxtl>re 

enthalpy  of  mixture 

enthalpy  of  water 

frozen  specific  heat  of  mixture 

specific  heat  of  vapor  at 


Assuming 


h  -  h 
m  w 


p,  m  p»v,f 

L  -  1000  Btu/lb 
*  5600  Btu/lb 


kf  «  2.  2  X  10'3  Btu/(ft)(sec)*K 


p,  w 


0.  24  Btu/lb 
62.4  lb /ft3 


and 


(H-12) 


(H-l  3) 


H-7 


we  readily  obtain 


t  =  4,  9  X  10'5(1  -  f  2/3} 
£ 


(K-  14) 


This  approximate  formula  leads  to  times  t  of  0.  032,  0.  038,  and  0.  049  msec 
for  80,  90,  and  :00%  evaporation.  The  application  of  these  results  and  the 
results  in  Appendix  H-5-a  are  discussed  in  Volume  I,  Section  III-E-2-b. 

6 .  DROPLET  DYNAMICS 

The  equation  cf  motion  of  a  droplet  parallel  to  the  air  stream  is  taken  as 


-  m  (U  .  -  U ,  ) 

dt  air  drop 


1  r  ltd 

2  UDPair  4 


(Uair-Udrop>  (H-'5: 


where  the  notation  is  standard.  If  the  water  jet  is  perpendicular  to  the  flow, 
x  is  the  distance  along  a  streamline,  y  is  the  distance  perpendicular  to  the 
streamline,  and  U.  is  the  initial  velocity  of  the  droplets,  then 


dx 

dt 


=  U  . 

air 


1  + 


D  P 


drop 


(H- 16} 


U. 

i 


(H-l  7) 


For  U  .  =15  ktt/sec,  =  2,  p  .  /p,  =  4  X  10'^  d  =  lOp, 

air  D  air  drop  n 

velocity  of  150  ft/sec,  the  following  results  are  obtained 


and  a  lateral 


t,  psec 

x,  cm 

y,  cm 

20 

0.  55 

0.  1 

100 

9.  7 

0.  5 

200 

31.2 

1.0 

These  values  are  the  basis  for  the  conclusion  that  the  lateral  spreading  of  the 
insectant  is  not  severe  (Volume  I,  Section  III-E-2-b). 
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APPENDIX  I 


ELECTROPHILIC  SEEDING 


Communication  difficulties  created  by  the  plasmas  that  surround 
spacecraft  during  atmospheric  reentry  have  stimulated  considerable  research 
on  the  effect  of  electronegative  gas  injection  on  the  ionization  levels  in  hot 
gases.  The  experiments  of  Carswell  and  Cloutier  (Ref.  1-1)  indicate  that 
concentration-  of  SF^  on  the  order  of  the  initial  electron  concentration  are 
very  effective  in  quenching  argon  plasmas  produced  by  an  rf  discharge. 
Conversely,  the  experiments  of  Fuhi  (Ref.  1-2)  and  Hoffman  and  Westbrook 
(Ref.  1-3)  indicate  that  SF^  is  much  less  effective  in  argon  arc  jets  and  seeded 
nitrogen-oxygen  flame  plasmas.  Theoretical  results  i. '  this  area  have  been 
either  specialized  (Refs.  I and  1-4)  or  inconclusive  (Ref.  1-3).  It  will  be 
shown  here  that  equilibrium  ionization  levels  of  high -temperature  air  are  not 
significantly  affected  by  the  presence  of  electronegative  species  but  that 
nonequilibrium  levels  can  be  rapidly  driven  to  equilibrium  levels  by  the 
electron-attachment  and  charge-transfer  mechanism.  These  conclusions 
are  consistent  with  the  experimental  observations,  and  they  provide  the  basis 
for  evaluating  the  usefulness  of  these  mechanisms  in  Sections  III-E,  IH-F, 
and  III-G  of  Volume  I. 

1.  EQUILIBRIUM 

The  equilibrium  concentrations  of  the  various  species  in  seeded  air 
for  a  given  temperature  and  pressure  are  obtained  by  simultaneously  solving 
the  chemical  equilibrium  equations  associated  with  the  pertinent  reactions 
and  the  appropriate  conservation  equations.  At  the  temperatures  and  pressures 
of  interest,  the  dominant  species  in  air  are  N2»  N,  O^,  O,  NO,  NO+,  and  e~. 

If  X  represents  the  lect ronegative  component,  the  mixture  will  also  contain 
X,X*,  and  whatever  compounds  can  be  formed  by  X  with  N  and  O.  The 
equilibrium  electron  density  will  be  a  function  of  temperature,  pressure,  and 
initial  seed -gas  concentration.  The  exact  solution  of  this  problem  requires 


I-i 


detailed  knowledge  of  chemistry  and  use  of  an  electronic  computer;  the 
problem  is  considerably  simplified,  however,  if  only  the  reduction  in  electron 
density  as  a  function  of  seed-gas  concentration  is  required.  If  the  optimistic 
assumption  is  made  that  the  electronegative  species  interact  only  with 
electrons,  and  the  reasonaole  assumption  is  made  that  the  variations  In  the 
concentrations  of  1*  and  O  due  to  the  injection  of  the  electronegative  gas  are 
negligible,  we  Obtain 


n  l  NO+J 

TnTToT 


ne(N£  -  Ne) 


N, 


(1-1) 


ne  +  NE  =  [NO4] 


where  n°  is  the  initial  electron  concentration,  n  is  the  equilibrium  electron 
e  e  .  —  . 

concentration,  Np  is  the  initial  electronegative  species  concentration,  NE  is  the 

equilibrium  negative  ion  concentration,  Kj  is  the  chemical  equilibrium  constant 

for  the  ionization  reaction,  is  the  chemical  equilibrium  constant  for  the 

electron  attachment  reaction,  EA  is  the  electron  affinity  of  the  electronegative 

species,  and  the  other  symbols  are  standard.  These  equations  lead  to 


/M2  "«  ♦  ka  >  ka 

\7)  =NE  +  n.+KA  NE  +  KA 


(1-2) 


1-2 


This  expression  differs  from  the  analogous  equation  of  Ref.  1-3  because  here 
account  has  been  taken  of  Lc  Chatelier'a  principle  —  as  the  electronegative 
component  removes  electrons,  the  reaction  N  +  O  —  NO+  +  e*  will  produce 
more  electrons.  Known  electron  affinities  do  not  exceed  4  V,  therefore  the 
limit  of  the  attachment  mechanism  effectiveness  can  easily  be  calculated. 
Such  calculations  for  10%  concentrations  of  4-V  electronegative  species  in 
air  have  been  performed  for  temperatures  and  pressures  characteristic  of 
the  reentry  environment,  and  the  results  are  presented  in  Fig.  1-1.  The 
curve  labelled  p/pQ  =  10"^  compares  favorably  with  exact  calculations 
(Ref.  1-4)  for  5%  concentrations  of  F^,  thus  indicating  the  validity  of  the 
assumptions  described.  The  electronegative  concentrations  required  for  a 
reentry  system  were  calculated  from  Eqs.  (1-1)  and  (1-2),  and  the  results 
are  presented  in  Volume  I,  Section  III-F-2. 

2.  NONEQUILIBRIUM 

When  the  electron  density  is  higher  than  the  equilibrium  level,  the 
approach  to  equilibrium  in  the  presence  of  an  electronegative  component  is 
given  by 


dn 

"dt 


•  •  -  «.,n2 
R  e 


o.N-n 
A  E  e 


(I-?) 


where  is  the  electron-ion  recombination  rate  constant,  and  a ^  is  the 

electron  attachment  rate  constant.  As  pointed  out  by  Carswell  and  Cloutier, 

we  may  assume  that  the  rate  at  which  the  electronegative  component  transfers 

electrons  to  positive  ions  is  much  faster  than  the  attachment  rate,  thus  N— 

0  ** 

may  be  regarded  as  a  constant.  If  we  make  the  transformation  t  T  *° 

that,  in  the  absence  of  the  electronegative  component,  the  electron  density 
decays  to  one -half  its  initial  value  when  t  =  !,  we  obtain  the  solutions  given 
in  Fig.  1-2,  where  the  seeding  parameter  r  represents  the  ratio  o^Nj./<r^n®. 
The  experimental  curves  presented  by  Carswell  and  Cloutier  are  in  agreement 
with  the  theoretical  curves  for  a  seeded  argon  plasma  when  the  reasonable 
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itmurm.H 


Fig.  I-t.  Reduction  of  Equilibrium 
Electron  Concentrations  in  Seeded  Air 


Fig.  1-2.  Reduction  of  Non- 
equilibrium  Electron 
Concentrations  in 
Seeded  Ionized 
Gases 
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-9 

assignment*  o0,  a  .  ~  10  are  made  (Refs.  1-5  and  6).  (The  value  for  ot.  is 

K  A  -16  2  A 
based  on  an  electron  attachment  cross  section  of  10"  cm  .  )  The  dashed  line 

in  the  figure  represents  the  measured  curve  for  the  unseeded  plasma  —  the 
discrepancy  is  probably  due  to  the  fact  that  the  electrons  must  be  cooled  be¬ 
fore  they  can  recombine  with  positive  ion*  (the  reported  initial  temperature 
was  7  x  104#K). 

It  may  be<concluded  from  these  considerations  that  electrophilic  seeding 
is  not  effective  in  reducing  equilibrium  electron  density  levels  in  nonexpanding 
flows  but  that  it  can  be  effective  in  catalyzing  the  electron-ion  recombination 
process  in  expanding  flows  or  in  flows  which  have  been  cooled  by  injection  of 
foreign  material.  The  specific  application  of  these  conclusions  to  the  reentry 
communications  problem  have  been  discussed  in  Volume  I,  Section  IJ2-E. 

III-F,  and  HI-G. 
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APPENDIX  J 


QUASI -OPTICAL  AND  OPTICAL  SYSTEMS 


The  power  requirements  for  the  systems  described  in  Volume  I, 
Section  III  - 1  arq  derived  here. 

1.  MILLIMETER  SYSTEM 

The  power  PR  required  at  the  input  of  a  state-of-the-art  95-GHz 
receiver  is  given  by 


PR  =  (S/NKNF)kTAf 


(J-i) 


where 

S/N  =10  (signal-to-noise  power  ratio) 

NF  -  40  (receiver  noise  figure) 

-21 

kT  =  4  y  10  W/Hz  (noise  power  per  cycle) 

Af  =  300  kHz  (receiver  bandpass) 

We  readily  obtain  PR  =  4.  8  X  10”* ^W.  The  required  output  power  PT  for  an 
omnidirectional  system  is  given  by 


PT  =  4tr  (R2/Ar)  PR  (J-2) 

where  R  is  the  range,  and  Ap  is  the  effective  area  of  the  receiver  antenna. 

n  2 

For  a  range  of  100  miles  and  an  area  of  200  ft  we  obtain  P^,  =  9  mW. 

Allowing  for  miscellaneous  losses  of  10  dB;  we  estimate  a  requirement  of  90 
mW  for  an  omnidirectional  clear-weather  system.  This  estimate  is  discussed 
in  Volume  I,  Section  III  - 1  - 1 . 


J-l 


2. 


INJECTION  LASER  SYSTEM 

The  signal-to-noise  ratio  at  the  output  of  a  photomultiplier  is  given  by 


ZeAfdg  +  +  Ijj) 


(J-3) 


where  Ig  is  the  current  due  to  the  signal,  Ig  is  the  current  due  to  the  back¬ 
ground  radiation,  and  1^  is  the  dark  current  of  the  tube.  In  this  case,  the 
dark  current  is  negligible,  and  the  required  received  power  is  given  by 


PR  =  GWq)(S/N)  Afjl  +11+  (2qPB/hvAf){S/N)'JJ1/2}  (J-4) 


where  PR  represents  the  background  radiation  power,  hv  is  the  photon 
energy,  q  is  the  quantum  efficiency  of  the  photomultiplier,  and  Af  is  the 
bandpass  of  the  video  system. 

The  noise  power  is  given  by 


PB  *  +  lV  <V*  >  ARB 


(J-5) 


where 


‘r 

ar  = 


3  12  1 

spectral  density  of  daylight  noise  *10”  W  sr  m  A” 

12  1 

spectral  density  of  vehicle  radiation  *10W  «r*  m”  A’ 

2 

receiving  aperture  *1  m 


B  =  optical  filter  bandpass  *10A 
«r  =  receiving  field  of  view  *10  ^sr 


using  the  indicated  estimates,  we  obtain 
Pjj  ®  10"^\V  (daylight  noise) 

Py  as  2  x  10  ~^W  (vehicle  noise) 

Pb«3*  10-8W 

-  3  _  i  q  g 

Taking  q  *  3  x  10  ,  hv  =  2.  2x10  J,  and  Af  =  10  Hz  (corresponding 

*  -8 
to  a  pulse  length  of  1  0  nsec),  we  obtain  P^*s  3x10  W.  If  atmospheric 

losses  are  estimated  as  6  dB,  a  requirement  of  50  kW  pulses  for  an  omni¬ 
directional  transmitting  system  with  a  range  of  100  miles  is  indicated.  If 
the  transmitted  power  can  be  concentrated  in  a  beam  with  a  width  of  about 
1  0  deg,  the  peak  power  requirement  is  reduced  by  a  factor  of  about  500, 
indicating  a  requirement  for  100-W  pulses.  As  discussed  in  Volume  I, 

Section  ni-I-2,  this  estimate  is  a  reasonable  requirement  for  future  systems. 
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